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Abstract 

We formulate an effective medium (mean field) theory of a material consisting of randomly dis- 
tributed nodes connected by straight slender rods, hinged at the nodes. Defining novel wavelength- 
dependent effective elastic moduli, we calculate both the static moduli and the dispersion relations 
of ultrasonic longitudinal and transverse elastic waves. At finite wave vector k the waves are dis- 
persive, with phase and group velocities decreasing with increasing wave vector. These results are 
directly applicable to networks with empty pore space. They also describe the solid matrix in two- 
component (Biot) theories of fluid-filled porous media. We suggest the possibility of low density 
materials with higher ratios of stiffness and strength to density than those of foams, aerogels or 
trabecular bone. 

PACS numbers: 43.20.Hq,43.20.Jr,43.35.Bf,43.80.Cs,87.10.Pq 



Electronic address: katz@wuphys.wustl.edu 



I. INTRODUCTION 

Many materials of biomedical or technological interest [l| contain a network of slender 
elastic rods with hinged connections at irregularly or randomly distributed nodes. In some 
cases the material may consist only of the network, while in others, such as trabecular 
bone, the space between the rods may be filled with fluid. Calculation of the elastic and 
acoustic properties of such a random network requires a statistical model. There have been 
extensive studies of fluid- filled networks ( (2- 7 ] ) but these have generally concentrated on 
the interaction between the matrix and the pore fluid, while treating the matrix as a simple 
homogenous material. 

The mechanical properties of networks of elastic rods are sensitive functions of their 
geometry. For example, simple cubic lattices have bulk moduli that scale as the volumetric 
filling factor J 7 of the rods but have zero stiffness to shear along the lattice planes, while 
in FCC and BCC lattices all moduli scale oc J 7 . The moduli of open-celled foams scale 
oc T 2 while those of closed-cell foams scale oc J 73 [lj]. In general, if a deformation requires 
extension or compression of rods the corresponding modulus scales oc J 7 , but if it is resisted 
only by rod flexure or by the resistance of the nodes to changes in the rod angles the moduli 
scale as a higher power of J 7 . Foams are soft because they can deform by bending their thin 
members without extension, changing angles at nodes or cell junctions with little resistance. 

Here we develop an elementary effective medium model of the elastic properties of a 
statistically isotropic matrix consisting of thin elastic rods joined by hinges at randomly dis- 
tributed nodes. This is surely an oversimplified description of a real material, as any analytic 
model of a random material must be, but may provide a useful guide to and parametriza- 
tion of its properties. Statistical modeling is necessary because the full three-dimensional 
structure of a solid consisting of irregularly located and connected nodes is unlikely to be 
known quantitatively and likely varies from realization to realization. In a biomaterial like 
cancellous bone it may vary from individual to individual, and within a single bone in a 
single individual. Even in statistically homogeneous foams it depends on details of their 
preparation. 

By connecting nodes of mean coordination number C > D, where D is the dimensionality, 
with straight rods that terminate on the nodes, we construct a model in which (except for 
pathological cases like a simple cubic lattice) strain implies changes in the distances between 
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nodes that are first order in the strain, and hence extensions of the rods of that order. The 
resulting structure is stiffer than open or closed cell foams such as expanded polymers, 
aerogels, hydrogels, sponges and spongy (trabecular or cancellous) bone in which rods or 
cell edges join at points along their lengths jl|. In our model stiffness is derived from the 
resistance of the rods to extension, rather than to bending; the former is much greater 



for slender rods. In contrast to structures near the isostaticity limit |8HlO|, the nodes are 
overconstrained. The great stiffness of structures like those we model may have practical 
utility. 

Within this model we calculate the dispersion of longitudinal and transverse elastic waves 
at finite wave vectors, where heterogeneity is significant. In this non-dissipative and non- 
local model dispersion is a consequence and measure of the inhomogeneous spatial structure 
(nonlocality invalidates one of the assumptions of the Kramers- Kronig relations). The dis- 
persion relations, functions of a single scalar variable (the magnitude of the wave- vector) , 
are a compact description of a complex three-dimensional structure. 

II. THE MODEL 

The model consists of nodes randomly, but statistically uniformly, distributed in space 
with a mean density n. In contrast to the related problems of the description of vibrations 
in a solid in which massive atoms are connected by massless elastic rods or springs [111] , or 
a granular material in which grains are pressed into frictional contact Q, in our model the 
nodes are massless and serve only as geometric constraints on the rods connecting them. 
Each node is connected to an average of C (the coordination number) nearest neighbor nodes 
by uniform and identical (except for length) elastic rods that are hinged at the nodes. We 
expect the discrete random structure to approach the continuum (effective medium) model 
as C — > oo. We present results for C = 4 and C = 12 to show the dependence of the model 
on C. Even C = 4 is well above the matrix's isostaticity threshold 8-loj], and is sufficient 
to assure its stiffness against shear. 

The distribution of distances ii to the z-th nearest neighbor of a node is obtained from 
an Erlang distribution (Gamma distribution for integer index) for the volume Vi enclosed 
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250 nodes/volume; C=12 








FIG. 1: Rods connecting randomly distributed nodes. One realization of randomly distributed 
nodes with mean density 250 within a unit cube is shown. The closest nodal neighbors are connected 
by rods, for a mean coordination number C = 12; some of these are outside the cube. The network 
is characterized by a length n -1 / 3 = 0.159 and is not fractal. The complexity of spatial structure, 
including voids and tongues of multiply connected nodes, is evident, but difficult to quantify except 
by its mechanical properties. 

by a sphere of radius ti containing its % nearest neighbors: 
The probability distribution of ii is 

P i (£ t ) = 4n^P v (Vi), (2) 

where £ { = (3^/4^) V 3 . 

Fig. 1 shows one realization of a network of rods in a unit cube connecting randomly 
distributed nodes. The nodes have a mean density n = 250 and mean coordination number 
C = 12. 
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In a microscopic realization of randomly distributed nodes, such as that shown in Fig. 1, 
it is not possible to enforce the condition that all nodes have the same coordination number; 
if node f is the C-th nearest neighbor of node j, to which it is connected, node j may (for 
example) be the C + 1-st nearest neighbor of node j'. Such complications are ignored, by 
definition, in an effective medium theory, in which C may be regarded as a mean coordination 
number. 

We consider the dispersion relations of elastic waves in this medium of thin straight 
massive elastic rods connecting massless nodes, with the rods completely hinged (free to 
change direction and to rotate about their axes) at the nodes. The effective medium model 
describes the motion of the nodes by a sinusoidal plane wave. The latter approximation is 
valid in the limit kn^ 1 ^ 3 — > 0, but breaks down for short waves with kn~ 1 ^ > it, for which 
scattering and localized modes are also important [13I [ijj]. In the limit of small amplitude 
(an assumption made throughout) compressive loads are below the rods' Euler buckling 
thresholds. 

The fundamental equations for a small amplitude plane longitudinal wave propagating 
in the x direction in an infinite medium are 



dv x 1 da xx 



dt p dx 
9« TT dv„ 



dt dx 



X 

1 



(3) 
(4) 



where a xx is the xx component of the stress tensor, u xx <C 1 is the xx component of the 
strain tensor describing the elastic wave, v x is the x component of material velocity and po 
is the mean density (including void space) of the matrix. In an effective medium model we 
consider a xx , u xx and v x to be continuous functions of space, as they would be in a continuous 
medium. The wave variables v x /cq and u xx /Bi{k) are of the same order of smallness as u xx , 
where Co is the thin fiber longitudinal sound speed of the rods and Bi{k) is the effective 
longitudinal modulus (defined below) of the matrix. 

These equations are closed with a constitutive relation 

o~xx = Bi(k)u xx , (5) 

where the effective longitudinal modulus Bi(k) is independent of space in an effective medium 
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theory, but depends on the mechanical properties of the rods and is a nontrivial function 
of the their spatial structure and of the wave vector k (a scalar in an isotropic effective 
medium). By taking Bi(k) to be real the model is lossless by assumption. A dissipative 
model could be defined by taking Bi(k) to be complex or, equivalently, by adding a term 
B[(k)u xx to Eq. 

The modulus is defined by the relation 

£ el = ^Bi(k)u 2 xx , (6) 

where the elastic energy density £ e i is obtained by calculating the elastic energy added to 
the rods by the strain field u xx (x,t). 
The resulting wave equation 

d 2 v x Bi{k) d 2 v x _ Q ^ 



has solutions 



with the dispersion relation 



dt 2 p dx 2 



v x0 exp i{kx ± ait), (8) 



For infinitesimal displacements of a hinged straight rod's endpoints it remains straight 
and below its buckling limit, so that its elastic energy is that of stretching. For slender rods, 
bending and torsional energy may be ignored even if the nodes are not hinged. Taking a 
cross-section A, Young's modulus E and equilibrium length £, stretched by an amount (not 
positive-definite) A£ satisfying \A£\ £, the elastic energy is 

ft- = ^. do) 

We describe a longitudinal elastic wave of amplitude a propagating in an infinite medium 
in the x direction with wave vector kx by a displacement field 

u x (x) = asm {kx + £), (11) 
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corresponding to a strain field 

Uxx = ok cos (kx + £)■ (12) 

The difference in ^-displacements of the ends of a rod, one end of which defines x = and 
the other is at x = £cos8 + 5u x , where 9 is the angle between the unstrained rod and the x 
axis, is 

5u x = a[sin (Hcos# + () — sin£]; (13) 

(^My = 5u z = 0. The choice of one end of the rod as defining x = is equivalent to a choice 
of phase £; either may be chosen freely if, in the computation of total elastic energy and 
modulus, an average is taken over the other. 

The rod undergoes a change in length (to first order in a) 

A£ = \/l 2 sin 6 2 + [I cos 6 - 5u x ) 2 - i (14) 

w cos 95u x (15) 

~ a cos #[sin (hi cos # + C) ~ sm C]- (16) 

Its elastic energy is 

AE AE 

£rod = -^j- cos 2 65u 2 x = "^jr® 2 cos2 $ [ s i n cos ^ + C) ~~ sm C] • (17) 

The mean elastic energy per volume is 

C AEa 2 I cos 2 6 r . „ „ >N . >n o\ 1 „ s / ? \ 

£d = x - r~ [sm ( + c) _ smC] / = 2 /( } { ' ( } 

where the factor of allows for the presence of two rods per node of coordination number 
C. 

For the sinusoidal wave ( TT2l) 

(u 2 xx ) = \a 2 k 2 (19) 
B t (k) = nC ^(^p?- [sin (hi cos 9 + C)-sm(} 2 \ . (20) 



and 

k 2 
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The dispersion relation (Q becomes 



- Jyl^/^tsm^cose + O-sinC] 2 ) , (21) 

\Z)Pm \ <- / 

where po = n^A(£)p m and p m is the material density of the rods. The coordination number 
enters only through the average over £, which is weakly dependent on C at finite kn^ 1 ^ 3 . 
The long wavelength limit is: 

JmftW -» f S, (22) 
independent of C, where the volumetric filling factor 

T= = \nCA{i). (23) 




The dispersion relation in this limit is 

^ = vf ' (24) 

where tVod = y/E/ p m is the phase and group velocity of the longitudinal wave in an indi- 
vidual rod in the (non-dispersive) thin- rod limit. 

The derivation has assumed J 7 <C 1, a condition implicit in the description of the matrix 
as a network of long slender straight rods with £ A 1 / 2 . If this condition is not met, then 
the assumption that the rods do not intersect between the pre-determined nodes is invalid. 



III. TRANSVERSE DISPERSION RELATION 



The dispersion relation of transverse waves is found from an analogous calculation. Defin- 
ing y as the polarization direction, we replace u x by u y , v x by v y , u xx by u xy and o xx by 
a xy . In Eqs. [HI HH dU 1201 and I2T1 outside the brackets sin9 is replaced by cos6> and cos# 
by sin# cos0, where is the azimuthal angle of the rod, taking = in the x-y plane. 
Averaging over introduces additional factors of \. The results Eqs. [22} EH for the static 
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TABLE I: Elastic constants (from Eq. [TKj) and wave speeds for propagation along crystal axes for 
cubic lattices of rods. 

limit bpcomp 



KmBtik) ->• ^-TE (25) 
fe^o 15 



W^vi^' (26) 



IV. RESULTS 



The dispersion relations for the longitudinal and transverse modes are shown in Fig. 2 
for C = 4 and C = 12. The frequency dependences of the corresponding phase and group 
wave velocities v p h and 

Vgr «3X6 shown in Fie;. 3. 
From the preceding results for the longitudinal and transverse wave speeds the static 
(u, k — >■ 0) Young's modulus E matrix and Poisson's ratio v m atrix of the bulk matrix may be 



obtained 



15|: 



E m atrix — n ^~E (27) 


^matrix J- (28) 

The model predicts the ratio v t r/vi ong = 1/\M of long wavelength (A ^> rir 1 / 3 ) wave speeds, 
independent of the properties of the rods. 

For comparison, we present results for the elastic constants of cubic lattices of rods and 
for the propagation speeds along the crystal axes of long- wavelength waves in Table HI 

The acoustic properties of fluid-filled porous media, such as aquifers, petroleum reservoirs 
and cancellous (trabecular) bone, are described by models originally developed by Biot jj- 
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Dispersion Relation 




1 



FIG. 2: Dispersion relations u(k) of longitudinal and transverse waves. v ro d = \jEj p m is the 
slender rod longitudinal wave propagation speed of the individual rods making up the matrix. The 
effective medium theory is expected to break down for fcn -1 / 3 > ir where scattering is strong. At 
finite /cn™ 1 / 3 these dispersion relations depend on the coordination number, as shown for C = 4, 12. 
Dispersion is greater for larger C because of the presence of longer rods with larger values of k£; 
for C = 4, (^n 1 / 3 = 778 while for q = 12) <^) n V3 = 10 85. 

Q|. These are also effective medium models that describe the wave modes of the fluid 
filled medium as if it were homogeneous and continuous, and do not attempt to include 
the dispersion at finite k, other than that which is a consequence (by the Kramers-Kronig 



relations 



16Nl8|) of dissipation. The properties of the empty matrix, such as those we have 



modeled, are required to calculate the properties of the Biot modes. 

We have found elastic moduli (Eq. 123 and Table [[]) that are proportional to the density 
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FIG. 3: Phase and group velocities of longitudinal and transverse waves as functions of fre- 
quency. C is indicated in subscripts. The effective medium theory is expected to break down 
for (jj/{p}' 3 v roc i) > 1. At finite /era -1 / 3 these velocities depend on the coordination number. 



or filling factor. Our model shows how stiff low density structures can be, if designed or 
evolved to maximize stiffness. Our model describes a random structure; similar stiffness 
has been found for ordered structures This is in contrast to aerogels [2( 
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25j that 



typically show E oc with 2 < m < 4. These values are characteristic of closed-cell foams 
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found m f« 1.5. The structures of 



2l|-|25| 



for which m — 3 although one study of aerogels 
aerogels are complex and sensitive to their preparation 

A study of water-filled trabecular bone \z\ found results that can be fitted by a value 
for the matrix of m ^ 3.7. These properties are very different from those of our model, 
and mechanical stiffness cannot have driven the evolution of trabecular bone. The remark- 



ably low modulus of the trabecular matrix explains the observation [27| that at low matrix 
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filling factor the fluid properties largely determine the elastic wave speed and attenuation 



fluid-filled trabecular bone have 
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in trabecular bone. This may explain why Biot models o 
required phenomeno logical fitting of their parameters 28- 

In our effective medium model dispersion is a consequence of microstructure and is present 
even without dissipation. The apparent contradiction with the Kramers- Kronig relations as 
adopted in [16 



crystals jll| , by noting that the effective medium model, like the phonon models, is nonlocal: 
the force on a mass is determined by the instantaneous positions of other masses. This 
violates the Kramers-Kronig assumption of locality, that response at one point is determined 
only by the history of fields at that point 33] . The mechanical compliance, analogous to the 
dielectric permittivity, is here a function of both u and k. Our model makes first-principles 
quantitative predictions of the dispersion and elastic wave speeds of vacuum-filled matrices 
that may be tested experimentally. 



18j is resolved, as it is in textbook point mass and spring models of phonons in 
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